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■^ ' Abstract 
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^S^ • High energy v^ and z>e propagation through the Earth has been performed using the 

'j^. Monte Carlo technique. We focused our attention on the effect of neutral current 

Mh. deep inelastic interactions ui N — > ui X compared to that of charged current ones 

,^ . ui N — > I X. We have found that NCs do not produce any significant effect with 

respect to the case in which only CCs are considered. Therefore we conclude that 
NC interactions can be neglected without considerable loss of accuracy. When com- 
puting upward-going neutrino fluxes a simple formula describing the transmission 
5^ I probability, that depends on the neutrino direction and energy and the CC cross 

section, can be used to account for the Earth shadowing effect. 
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1 Introduction 



Numerous works were published during tlie last decade on high energy neu- 
trino propagation through the Earth (see, e.g. [1-11], though the whole list of 
publications dedicated to this item is much longer). This is an important topic 
since, in order to correctly estimate the event rates due to neutrinos crossing 
the Earth that interact close or inside underwater neutrino telescopes (UNTs) 
[12-17], the Earth shadowing effect needs to be accounted for. 

The most complex case is that of Vr- Besides suffering neutral current (NC) 
interactions, it undergoes a regeneration chain of processes in the Earth due 
to charged current (CC) interactions and consequent fast r-lepton decay, pro- 
ducing a Pr of lower energy than the initial one: Vr — ^ r ^ z/^ . . .. More- 
over, secondary z/g's and z/^'s are generated via r+ -^ Ugi/re'^, r^ -^ UgiyrC', 
r+ -^ Ufj^Ur fJ''^ and r^ — ^ iy^Urfi", decays. Hence, a Ur with initial energy at 
Earth of E^^ leaves it with an energy El^ < E^^ , being accompanied by some 
amount of secondary ly^^eS [5,8]. 

In contrast to z^^'s, z^g's and t'^'s can only loose energy undergoing a NC 
interaction ui N — > ui X or be absorbed due to a CC interaction z// A^ — > 
IX ^ . Generally, to account properly these two processes one has to solve 
complex transport equations or apply Monte Carlo (MC) simulation. However, 
if NC interactions are neglected, the neutrino flux emerging from the Earth 
^i{Ey^^,0) is given by the following simple formula which depends only on 
the neutrino energy Ey^^ ^ and nadir angle 9: 

$i(£;,^,„^) = 5(i?.„„^) X %{E,,,^,e). (1) 



Here ^q{Eij ^,9) is the incident flux and S{Ei, ^,6) is the screening factor (or 
transmission probability) which is deflned as [1]: 



S{E,^^^,9) = I -paUE.,,^,0) = exp ]^-NAXacciE,^Jy<l piOJ) dl 



,(2) 



where pabs{Ey^^,9) is the absorption probability, acc{Eyf,.e) is the total CC 
cross section, Na is the Avogadro number, p{9, 1) is the Earth density proflle 
for a given direction 9 and / is the distance in the Earth on which the integral 
is performed to calculate the column depth. 

In this short paper we describe the results of the MC simulations for z/„ and 



^ We do not consider W~ production through the resonant process i>e e~ -^ W~ -^ 
X at £^pg«6.3PeV (the Glashow resonance [18]) which occurs only to u^, in a very 
narrow energy range. 



Ve propagation through the Earth, focusing on the role of NC interactions and 
analyzing the accuracy of the formula in Eq. (1). 



2 Method and Results 



The simulation we have used is fully described in [8]. Shortly, we remind that 
deep inelastic cross sections are calculated using CTEQ3-DIS parton distri- 
bution functions [19] taken from PDFLIB library [20]. We have assumed the 
Earth composition made by standard rock {A=22, Z=ll) of variable density 
with the Earth density profile taken from [21] (Fig. 1). 
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Fig. 1. The Earth density profile [21] (upper panel) and the Earth column depth 
vs nadir angle (lower panel) as it is used in our calculations. 



2.1 Results on monoenergetic neutrino beams 



We have simulated 1640 monoenergetic neutrino beams each of 5x10^ events 
(z^e, i/e, y^ and v^, 41 neutrino energies in the range 10 GeV < E^ < 10^ GeV 
and 10 nadir angles in the range 0.1<cos^<1.0) and we have saved informa- 
tion on neutrinos which survive after propagation. 

The CC cross section crcci that provides the probability of muon and elec- 
tron neutrino absorption, grows with energy and consequently the interaction 
length L^*^ = (A^^ x cfcc)~^ (in water equivalent units) decreases. At some 
'transparency' energy L^'^ becomes equal to the Earth column depth seen by 
the neutrino for a given nadir angle 9 and so, at E^^^ > Etransp{d) the Earth 
begins to consistently absorb I'e./i's. Etranspi^) grows with the nadir angle. It 
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Fig. 2. Left upper panel: CC interaction length vs E^^^^. The three horizontal lines 
are the Earth column depth for 3 different values of the nadir angle {cos 6 = 0.1, 
0.5, 1.0). Left lower panel: transparency energy vs cos0. The transparency energy 
is defined as the energy at which the neutrino interaction length equals the Earth 
column depth. A 'knee' at cos^ ~ 0.8 corresponds to an 'ankle' at the curve for 
the Earth column depth vs nadir angle (Fig. 1) which originates from the fact that 
the Earth density changes sharply at the distance ~3500km from the center. Right 
panel: absorption probability for f^^e's {u^^es) in the Earth. In the simulation of 
propagation through the Earth only CC interactions were taken into account. In all 
panels solid lines are for neutrinos and dashed ones for anti-neutrinos. 



0.8 



tO.6 



0.4 



0.2 




0.12 



0.1 



£ 0.06 

Q. 



0.04 



0.02 



3 4 5 6 7 

log,o(Ev. GeV) 



O cos(e): 

cos(e): 
cos(e): 
O cos(e): 




3 4 5 6 7 

log,„(E^, GeV) 



Fig. 3. Left panel: neutrino absorption probability in the Earth taking into account 
(dashed lines) or not (solid lines) NC interactions. Right panel: probability for fe,^'s 
(Pc,^'s) to pass through the Earth with no CC interactions but with at least one 
NC interaction (solid lines are for neutrinos and dashed ones for anti- neutrinos) . In 
both plots results are given for four nadir angles {cosO = 0.1, 0.4, 0.7, 1.0). 



is about 50-80 TeV at directions close to the vertical {9 = 0°), it becomes 
about 60 PeV at cos(^) = 0.1 (^=84°) and increases to larger values for more 
horizontal directions since the matter to be transversed decreases (Fig. 2, left 
panels). For very horizontal directions the composition of the Cherenkov me- 
dia where the array is located or also the details of the coast for a specific 
underwater detector should be known to calculate exactly the transparency 
energy. Anyway, its value will be very large. Correspondingly, the probability 
that a neutrino undergoes a CC interaction and is absorbed grows with E^, ^. 
In the right panel of Fig. 2 the absorption probability pabsiE^^,^^, 9), computed 
as the ratio between the number of absorbed and propagated neutrinos for a 
given energy and nadir angle, is shown in the case where only CC interactions 
are simulated. 

In the considered energy range, (Ttvc is 2.3-3.2 times smaller compared to ace 
and hence, the probability that a CC interaction occurs is higher than that of 
a NC. When fe,/i's interact by a NC, the neutrino energy decreases on average 
by a factor 1 — (?/) = 0.5-0.8 (depending on energy), since the rest of the energy 
{y)Ei, ^ is taken by the hadronic shower. Hence, the probability that another 
CC interaction occurs decreases. In other words, the probability pabs{E^^^^,9) 
that the neutrino will be absorbed becomes lower. But quantitatively the effect 
is very small as it can be seen in the left panel of Fig. 3 where pabs{Ep^^, 9) 
for z/e,^'s is plotted for propagation in the Earth with- and without accounting 
for NC interactions (results for z/e,/i's are qualitatively the same). Comparing 
the NCH-CC case to the CC only case, we observe only a slight decrease of the 
absorption probability. The right panel of Fig. 3 shows the probability p]^c that 
t'e.^'s survive after crossing the Earth with no CC interactions but at least one 
NC interaction. The maximum value oipNC for different neutrino nadir angles 
is about 0.13 for energies around the transparency energy Etransp{9). At lower 
energies it falls down due to the cross section aj\fc almost linear decrease with 
energy. At higher energies it is suppressed due to the increase of ace since, 
even after loosing energy in NC interactions, the neutrino energy remains high 
enough to undergo another CC interaction that makes the neutrino disappear. 
In any case, pnc is remarkably lower compared to pcc and, therefore, it does 
not affect significantly the total Pabs{E^^^, 9), which turns out to be very close 
to Pcc- 



2.2 Results on various energy spectra 



We have investigated the effect of muon neutrino propagation through the 
Earth ^ for various models of neutrino fluxes incident on the Earth when only 



^ All the results and conclusions relevant to v^^s are also valid for z/g's since differ- 



ential cross sections are practically the same. 
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Fig. 4. Upper row: muon neutrino spectra before (solid line) and after propagation 
through the Earth in the vertical direction with- (dotted lines) and without (dashed 
lines) accounting for NC interactions. Left: AGN Protheroe spectrum [22]; middle: 
MPR upper bound on diffuse neutrino flux [23]; right: atmospheric muon neutrino 
spectrum in [24] including both neutrinos from pion and kaon decay and prompt neu- 
trinos from charmed meson decay (calculated in the Recombination Quark-Parton 
Model frame). We considered the same amount of neutrinos and anti- neutrinos 
{4>u^ : 0p^ = 1 : 1). Lower row: the screening factor S{Ey^) = ^i{Ey^)/^Q{Ey^) 
computed with fluxes presented in upper plots. Open circles: accounting for CC 
interactions only; full circles: MC including both NC and CC interactions. 

CCs are simulated and when also NCs interactions are accounted for. 

We considered the spectrum predicted by Protheroe [22] for diffuse fluxes 
of neutrinos produced by active galactic nuclei (AGNs), a flux equal to the 
Mannheim, Protheroe and Rachen (MPR) upper bound [23] and the atmo- 
spheric muon neutrino spectrum from [24]. Results are shown in Fig. 4 only 
for the vertical direction, that is the direction where interactions are most 
probable given the largest column depth. Results on screening factors (shown 
in the three lower panels) are in agreement with the results reported in [3], 
but, as we will see below, large values of screening factors (S" > 1) cannot be 
interpreted as a significant enhancement of the emerging neutrino spectrum 
due to NC interactions (as it is concluded in [3]). As a matter of fact, this 
effect is due to the spectral shape. Considering upper panels in Fig. 4 where 
emerging neutrino spectra are given one can conclude that effects of NC in- 
teractions can be summarized in two points. Firstly, NCs result in a slight 
shift of the mean energy of the spectra emerging after propagation to lower 



energies compared to the spectra obtained simulating only CCs (the spectral 
shapes are almost the same). This energy shift decreases with increasing nadir 
angles 6. For the AGN spectrum [22] is equal to 22%, 9% and 3% for ^=0°, 
^=60° and ^=80°, respectively. The shift is due to the decrease of neutrino 
energy in NCs since a fraction of the parent neutrino energy is taken away 
by the hadronic shower, hence the resulting neutrino has lower energy than 
its parent. Secondly, some enhancement of emerging spectra is observed when 
NC interactions are accounted for with respect to the case in which the simu- 
lation is performed without NCs. The most remarkable enhancement is found 
for the AGN spectrum in Ref. [22] for which the ratios of the integral of the 
emerging spectra with and without including NC interactions are equal to 
1.13, 1.02 and 1.003 for ^=0°, ^=60° and ^=80°, respectively. Using the val- 
ues for ^=60° (cos^ =0.5) and this AGN flux to estimate the corresponding 
effect for the whole lower hemisphere (given that dependencies are more or 
less linear ones) we obtain a 9% energy shift in the spectral maximum and 
a 2% flux enhancement due to NC interactions. Thus, when NCs are taken 
into account, neutrino spectra are shifted to lower energies by an amount that 
hardly will be appreciated by UNTs for which the typical error on neutrino 
energy measurement is never less than a factor of 2 (see, e.g., [25]). Also the 
enhancement of neutrino fluxes (and consequently of counting rates) due to 
NCs is small compared to experimental and simulation uncertainties in UNTs 
[12-17]. 

The noticeable increase (up to S{Ey^) ^ 1.35) of the screening factor in the 
left bottom panel of Fig. 4 for the considered AGN spectrum [22] in the TeV 
range due to NC interactions is in a good agreement with [3]. But we believe 
that S{Eij ) is not the proper variable to estimate the effective enhancement 
of the emerging spectrum after propagation, and hence of event rates. The 
effect of NCs can be understood better if we compare the neutrino emerging 
spectra obtained with and without accounting for NCs. As shown in Fig. 4 
(upper panels) these spectra are very similar. Hence, the enhancement in terms 
oi S{Ey ) is the result only of the slight energy shift of the emerging spec- 
trum after propagation in the Earth relative to the initial spectrum at the 
Earth surface. The increase on the screening factor is larger for spectra with 
a maximum and steeply falling down for decreasing energies below the en- 
ergy corresponding to the maximum (see Protheroe AGN spectrum in Fig. 4). 
The effect is much smaller for flat spectra or power law spectra with negative 
spectral index, such is the case for MPR and for atmospheric neutrino spectra. 



3 Conclusions 



We have shown that the effect of NC interactions of ultra and extremely high 
energy muon and electron neutrinos propagating through the Earth on the 



neutrino flux emerging from tlie Eartfi does not play a significant role. Conse- 
quently event rates in UNTs, such as AMANDA [12], ANTARES [13], Baikal 
[14], IceCube [15], NEMO [16], NESTOR [17], are not much affected by NC 
effects during propagation in the Earth. Their effect can be neglected if the 
required accuracy of calculations is lower than ^10% in terms of energy mea- 
surement and few percents in terms of absolute event rates. Calculations of 
counting rates can be much faster if the simple formula in Eq. (1) is used 
compared to a full propagation of neutrinos through the Earth with MC sim- 
ulations or semi-analytical methods. 
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